Biopolymer Sequencing By Hybridization of probes to form ternary complexes and variable range alignment

ABSTRACT

Methods for sequencing a biopolymer by forming local ternary complexes along the length of the double-stranded biopolymer target molecule using one or more probes and obtaining information about the location of the probe(s) using a detector. These methods offer particular advantage when implemented with nanopore (including micropore) detection systems.

CROSS REFERENCE TO RELATED APPLICATIONS

This application is a continuation of U.S. Non-provisional patentapplication Ser. No. 12/243,451 filed on Oct. 1, 2008 now U.S. Pat. No.8,278,047, which claims the benefit of priority under 35 U.S.C. §119(e)to U.S. Provisional Ser. Nos. 60/976,714, filed on Oct. 1, 2007, and60/976,739, filed Oct. 1, 2007; the entirety of each application isincorporated herein by reference.

This application also incorporates by reference U.S. patent applicationSer. No. 11/538,189 entitled “HYBRIDIZATION ASSISTED NANOPORESEQUENCING”, which was filed Oct. 3, 2006, and published as U.S. PatentPublication No. 2007/0190542. This application also incorporates byreference U.S. Patent Application Ser. No. 61/093,885 entitled, “USE OFLATERALLY DISPLACED TRANSVERSE NANOSCALE ELECTRODES FOR VOLTAGE SENSINGOF BIOMOLECULES AND OTHER ANALYTES IN FLUIDIC CHANNELS,” by Ling et al.,which was filed Sep. 3, 2008.

TECHNICAL FIELD

The present invention relates generally to biopolymer sequencing. Moreparticularly, in certain embodiments, the invention relates to thedetermination of biopolymer sequence information using sequence-specificprobes capable of binding to double-stranded biopolymers. For example,positional information of one or more binding locations of the probesalong the target biopolymer molecule is determined and a sequence of thebiopolymer is derived based on the positional information.

BACKGROUND

There is great interest in identifying the composition and sequence ofvarious biomolecules, such as human DNA, with accuracy and specificity.Sequencing technology, however, is time consuming and expensive todevelop and implement. For example, sequencing the DNA of a singleindividual for the Human Genome Project required over $3 billion offunding.

It is estimated that each person's DNA varies from one another byapproximately 1 base in 1000. Knowledge of such genetic variations amonghuman populations may allow the scientific community to identify genetictrends that are related to various medical predispositions, conditions,or diseases, and may lead to the realization of truly personalizedmedicine where treatments are customized for a given individual based onthat individual's DNA. A reduction in the time and cost of DNAsequencing is needed to develop such knowledge and to tailor medicaldiagnostics and treatments based on the genetic makeup of individualpatients.

One particular obstacle inherent in known methods is the inability toaccurately position repetitive sequences in DNA fragments. Furthermore,known methods cannot determine the length of tandem short repeats, whichare associated with several human genetic diseases.

One emerging sequencing technology employs nanopore or microporedevices. Nanopores are substantially cylindrical holes formed in amembrane or solid media, said holes having diameters that range fromabout 1 nm to about 200 nm. Some existing methods using nanopores haveattempted to detect single DNA bases as they move through a nanoporeunder a bias voltage. However, it is difficult to detect single DNAbases as each base passes through the nanopore. Furthermore, the use ofnanopores small enough to track single stranded DNA are unreliable anddifficult to form.

Other methods have attempted to use nanopores to detect DNAhybridization probes or oligonucleotides on a DNA molecule and torecover the DNA sequence information using the method ofSequencing-By-Hybridization (SBH). SBH is a two step procedure, whereinthe collection of all subsequences that make up a target sequence isfirst determined by detecting hybridization of sequence-specific probesor a pool of probes to the target sequence and then using an algorithmthat relies on the use of combinatorial methods to reconstruct the fullsequence of the target using the collection of subsequences. Most of theSBH methods have relied on standard DNA probes, termed k-mers (see e.g.,E. M. Southern. “DNA chips: analysing sequence by hybridization tooligonucleotide on a large scale” Trends in Genetics, 12(3), 110-115(1996)).

SBH procedures can also be used to attach a large set of single-strandedfragments or probes to a substrate to form a sequencing chip. When asolution of labeled, single-stranded target DNA fragments is exposed tothe chip, the target fragments hybridize with complementary sequences onthe chip. The hybridized fragments can be identified using a radiometricor optical detector depending on the selected label. Each hybridizationprovides information about whether the fragment sequence is asubsequence of the target DNA. The target DNA can then be sequencedbased on which strings are and are not substrings of the targetsequence.

The efficiency of SBH methods is poor. For example, large probe arraysare required to sequence modest target lengths. Furthermore, theinformation regarding the binding position along the target sequence ofa given fragment with respect to other fragments is not generated usingthis experimental approach, and the number of times that a fragmentbinds a target is also undetermined. While SBH may be a useful forsequencing variants of known molecules, it is not useful for sequencingorganic biomolecules at high throughput and accuracy. Still further, thealgorithms that are used to reconstruct the target sequence from thehybridization data have not prove useful in practice because known SBHmethods do not return sufficient information to sequence long fragments.Thus, these limitations have prevented the adoption of SBH as a primarysequencing method. There is therefore a need for improved methods ofsequencing organic biomolecules such as DNA.

SUMMARY OF THE INVENTION

The present invention provides methods for sequencing a biopolymer byforming local ternary complexes along the length of the double-strandedbiopolymer target molecule using one or more probes and obtaininginformation about the location of the probe(s) using a detector. Thesemethods offer particular advantage when implemented with nanopore(including micropore) detection systems.

It is now appreciated that a key problem of analyzing single strandedbiopolymers through a nanopore is that single-stranded biopolymers tendto self-hybridize and fold into secondary or higher-order structures.This may prevent entry of the single stranded sequence into the nanoporeor stall its migration through the nanopore. Single-stranded DNA has ashort persistence length and a measure of the DNA polymer rigidity is ata local level. Thus, folding of single-stranded DNA may falsely providea signal similar to the presence of a binding probe.

Sequencing organic biopolymers comprising double-stranded DNA hasseveral advantages over the use of single-stranded DNA. Double-strandedDNA, in contrast to single stranded DNA, does not have a tendency tofold into secondary structures, is easier to manipulate prior to theintroduction into a nanopore channel, has a consistent transmissionspeed through nanopore channels, and has a longer persistence length. Inaddition, higher affinity probe molecules with greater sequenceselectivity may be used with double-stranded DNA, than withsingle-stranded DNA.

In one aspect, the invention features a method for determining a wholeor partial sequence of a target biopolymer. The method includes the stepof providing a double-stranded biopolymer target molecule. The targetmolecule (or fragment thereof) is contacted with a first probe havingspecificity for one or more first recognition sites of the targetmolecule, thereby forming one or more first local ternary complexesalong the double-stranded biopolymer target molecule (or fragmentthereof), the one or more first recognition sites each having a firstknown sequence. At least a subset of the one or more first local ternarycomplexes along the double-stranded biopolymer target molecule (orfragment thereof) is detected to determine a first probe map. The wholeor partial sequence of the target biopolymer may be determined using atleast the first probe map.

In certain embodiments, the detecting step includes passing thedouble-stranded biopolymer target molecule (or fragment thereof) througha nanopore and detecting an electrical signal indicative of one or morelocations of the one or more local ternary complexes along a length ofthe double-stranded biopolymer target molecule, thereby determining thefirst probe map.

The first probe map may include one or more of the following: (a)relative positional information of two or more first local ternarycomplexes along the double-stranded biopolymer target molecule (orfragment thereof); (b) absolute positional information of one or morefirst local ternary complexes along the double-stranded biopolymertarget molecule (or fragment thereof); and/or (c) error of positionalinformation of one or more first local ternary complexes along thedouble-stranded biopolymer target molecule (or fragment thereof).

In certain embodiments, the target molecule (or fragment thereof) iscontacted with a second probe having specificity for one or more secondrecognition sites of the target molecule, thereby forming one or moresecond local ternary complexes along the double-stranded biopolymertarget molecule (or fragment thereof), the one or more secondrecognition sites each having a second known sequence. The detectingstep may further include detecting at least a subset of the one or moresecond local ternary complexes along the double-stranded biopolymertarget molecule (or fragment thereof) to determine a second probe map.The determining step may include determining the whole or partialsequence of the target biopolymer using at least the first probe map andthe second probe map.

The first probe map and/or the second probe map may include one or moreof the following: (a) relative positional information of two or morefirst local ternary complexes along the double-stranded biopolymertarget molecule (or fragment thereof); (b) relative positionalinformation of two or more second local ternary complexes along thedouble-stranded biopolymer target molecule (or fragment thereof); (c)relative positional information of one or more first local ternarycomplexes and one or more second local ternary complexes along thedouble-stranded biopolymer target molecule (or fragment thereof); (d)absolute positional information of one or more first local ternarycomplexes along the double-stranded biopolymer target molecule (orfragment thereof); (e) absolute positional information of one or moresecond local ternary complexes along the double-stranded biopolymertarget molecule (or fragment thereof); (f) error of positionalinformation of one or more first local ternary complexes along thedouble-stranded biopolymer target molecule (or fragment thereof); and/or(g) error of positional information of one or more second local ternarycomplexes along the double-stranded biopolymer target molecule (orfragment thereof).

In certain embodiments, the determining step includes determining aspectrum map from at least the first probe map and the second probe map.Determining the spectrum map may include using a combinatorial sequencereconstruction algorithm. Determining the spectrum map may includedetermining a candidate sequence by ordering at least two probesequences using positional information and error of positionalinformation.

The first recognition site sequence and the second recognition sitesequence may overlap, and determining the spectrum map may includealigning the recognition site sequences. The first probe may be an n-merprobe, with n being an integer from 3 to 10. The first probe may includeone or more moieties, such as a polyamide, a nucleic acid, and/or apeptide nucleic acid. The first probe may include a gapped polyamidemolecule.

The contacting step may further include contacting the target molecule(or fragment thereof) with the first probe in the presence of at leastone of a recombinational protein or a recombinational enzyme so as tofacilitate hybridization of the first probe with the target molecule.

The providing step may include enzymatically synthesizing at least aportion of a second strand upon a single-stranded template. Thecontacting step may include binding the first probe to a minor groove ofthe target molecule (or fragment thereof).

The description of elements of the embodiments of other aspects of theinvention may be applied to this aspect of the invention as well.

In another aspect, the invention features a method for determining awhole or partial sequence of a target biopolymer. The method includesthe step of providing a double-stranded biopolymer target molecule. Thetarget molecule (or fragment thereof) is contacted with a plurality of nprobes each having specificity for one or more recognition sites of thetarget molecule, thereby forming local ternary complexes along thedouble-stranded biopolymer target molecule (or fragment thereof), eachof the recognition sites having a known sequence. At least a subset ofthe local ternary complexes along the double-stranded biopolymer targetmolecule (or fragment thereof) is detected to determine one or moreprobe maps. The whole or partial sequence of the target biopolymer isdetermined using at least the one or more probe maps.

In certain embodiments, the detecting step includes passing thedouble-stranded biopolymer target molecule (or fragment thereof) througha nanopore and detecting an electrical signal indicative of locations ofthe local ternary complexes along a length of the double-strandedbiopolymer target molecule, thereby determining the one or more probemaps. The detecting step may be repeated for each of the n probes.

The description of elements of the embodiments of other aspects of theinvention may be applied to this aspect of the invention as well.

In yet another aspect, the invention features a method for detecting thesequence of a target nucleic acid. The method includes contacting atarget molecule with a first probe to create at least one probe-targetcomplex at a recognition site of the target for which the first probehas a known specificity, while leaving uncomplexed, regions of thetarget for which the first probe is not specific. The target molecule iscontacted with a second probe to create at least one probe-targetcomplex at a recognition site of the target for which the second probehas a known specificity, while leaving uncomplexed, regions of thetarget for which the second probe is not specific. Complexed anduncomplexed regions of the target are detected and recorded to create afirst probe map of the first probe and a second probe map of the secondprobe, the first probe map and the second probe map incorporatinginformation on the relative position of the hybridization of the probes.A candidate sequence is determined by ordering at least two probesequences using positional information or a combination of overlappingsequences of the probe molecules and positional information.

One of more of the following features may be included. Detecting thefirst probe map and detecting the second probe map may be performedsequentially or concurrently.

At least one of the first probe and the second probe may include a tagso as to allow discrimination of bound first probe from bound secondprobe. The first probe and the second probe may be comprised of mixturesof nucleic acid sequences. Detection may include using a single moleculedetector, which may include a nanopore detector.

The description of elements of the embodiments of other aspects of theinvention may be applied to this aspect of the invention as well.

In another aspect, the invention features a method for detecting thesequence of a target nucleic acid. The method includes contacting atarget molecule with a first probe to create at least one probe-targetcomplex at a recognition site of the target for which the first probehas a known specificity, while leaving uncomplexed, regions of thetarget for which the first probe is not specific. The target molecule iscontacted with a second probe to create at least one probe-targetcomplex at a recognition site of the target for which the second probehas a known specificity, while leaving uncomplexed, regions of thetarget for which the second probe is not specific. Complexed anduncomplexed regions of the target are detected and recorded to create afirst probe map of the first probe and a second probe map of the secondprobe, the first probe map and the second probe map incorporatinginformation on the relative position of the hybridization of the probesand incorporating information on the error of the positional informationfor each probe. A candidate sequence is determined by ordering at leasttwo probe sequences using positional information and parameters relatingto the error in positional information or a combination of overlappingsequences of the probe molecules and positional information and error inpositional information.

One or more of the following features may be included. A parameterrelating to error in measured position of a probe may define a range ofpositions for the probe. The target sequence may be determined byoverlapping sequences of at least two probe molecules only when therange of one of the probes overlaps with the range of another probe. Therange may be selected so as to be proximal to a sequence position. Therange may be different for every probe. Choosing the range further mayinclude using a correlation between the size of the error and the sizeof the range. Detecting may include using a single molecule detector.The detector may include comprises a nanopore detector.

The description of elements of the embodiments of other aspects of theinvention may be applied to this aspect of the invention as well.

In yet another aspect, the invention features a method for determining asequence of a nucleic acid target. The method includes creating multiplealiquots of a target molecule solution. Each of a plurality of thealiquots is combined and stringently incubated with a differentcorresponding probe, with each probe-aliquot combination having aplurality of probe molecules having a known specificity for a givensubsequence of the target so as to create a probe-target complex foreach aliquot. The probes are predominantly complexed with theircorresponding subsequences for which they have specificity. Theprobe-target complexes are detected with a single molecule detectorcapable of distinguishing regions of the target that are complexed anduncomplexed with the probe so as to generate a probe map that estimates,to within an error, the positions of the multiple instances of thesubsequence within the target sequence. The probe maps generated fromthe multiple aliquots are preliminarily aligned, and a model sequencecreated that estimates the actual sequence of the target, the modelsequence having a growing end. A range of candidate sequences proximalto the growing end are selected from within the set of aligned probemaps. A probe-specific sequence is selected from within the range thathas an overlapping region of substantial homology or identity to asequence at the growing end of the model sequence, yet has at least onebase in a non-overlapping region beyond the growing end. At least onebase from the non-overlapping region is appended to the growing end.

One or more of the following features may be included. The range may beselected based on at least one estimated probe map positional error. Thesize of the range may be related to the error. The same range may beused to select sequences from each probe map. A different range may beused to select sequences from each probe map. The detector may be ananopore detector. Aligning the probe maps from the multiple aliquotsmay include using at least one of a reference probe, a referencesequence and a reference protein. The target may be fragmented prior tocreating the aliquots. Multiple fragment probe maps for a givensequence-specific probe may be aligned and assembled into a target probemap for that probe. The probes may be chosen so as to bind to allpossible sequences of DNA. Creating a model sequence may includechoosing a starting probe-specific sequence and from a probe map. Themultiple probe maps may be aligned in relation to at least one of eachother and the target sequence.

The description of elements of the embodiments of other aspects of theinvention may be applied to this aspect of the invention as well.

In yet another aspect, the invention features a method for detecting asequence of a target biopolymer comprising. The method includesproviding a double-stranded biopolymer target molecule. The targetmolecule is contacted with a first probe having a first probespecificity for recognition sites of the target molecule to form a firstplurality of local ternary complexes, the first probe having a firstpredicted recognition site sequence. A detector is used to determinepositional information of the first plurality of local ternarycomplexes.

One or more of the following features may be included. The positionalinformation may include a parameter related to a spatial distancebetween two local ternary complexes. The target molecule may becontacted with a second probe having a second probe specificity forrecognition sites of the target molecule to form a second plurality oflocal ternary complexes, the second probe having a second predictedrecognition site sequence, and a detector may be used to determinepositional information of the second plurality of local ternarycomplexes. The positional information of at least the first and secondplurality of local ternary complexes may be aligned to determine a DNAsequence of the target.

The description of elements of the embodiments of other aspects of theinvention may be applied to this aspect of the invention as well.

In yet another aspect, the invention features a method for detecting asequence of a target biopolymer. The method includes providing adouble-stranded biopolymer target molecule. The target molecule iscontacted with a first probe having a first probe specificity forrecognition sites of the target molecule to form a first plurality oflocal ternary complexes, the first probe having a first predictedrecognition site sequence. The target molecule is contacted with asecond probe having a second probe specificity for recognition sites ofthe target molecule to form a second plurality of local ternarycomplexes, with the second probe having a second predicted recognitionsite sequence. The formation of the local ternary complexes is detected.A sequence of the target biopolymer is determined by assembling at leastthe first and second predicted recognition site sequences.

One or more of the following features may be included. Detection mayinclude measuring the relative positions of the local ternary complexes,and/or creating at least one of a probe map or a spectrum map.Assembling may include overlapping the recognition site sequences. Oneof the first probe and the second probe may include a collection ofprobe molecules of substantially similar specificity. The probemolecules may include moieties such as polyamides, nucleic acids, and/orpeptide nucleic acids. Detection may include using a nanopore detector.Providing the double-stranded biopolymer target molecule may includeenzymatically synthesizing at least a portion of a second strand upon asingle stranded template.

The description of elements of the embodiments of other aspects of theinvention may be applied to this aspect of the invention as well.

In yet another aspect, the invention includes a library of probesadapted for determining a sequence of a biopolymer target, the probesincluding polyamide molecules. The library may also include gappedpolyamide molecules. All of the probes may be gapped polyamides.

The description of elements of the embodiments of other aspects of theinvention may be applied to this aspect of the invention as well.

In yet another aspect, the invention includes a method for determining asequence of an unknown sequence of a biopolymer target. The methodincludes: (i) contacting the target with a polyamide probe having aspecificity for a recognition site of the target; and (ii) detectingwhether the probe and target are bound.

One or more of the following features may be included. A predictedsequence of the recognition site may be recorded if the probe is bound.A probe map may be determined. The probe molecule may bind to the minorgroove of the target molecule. The description of elements of theembodiments of other aspects of the invention may be applied to thisaspect of the invention as well.

In yet another aspect, the invention includes a method for determiningsequence information of a nucleic acid having a sequence. The methodincludes providing a double-stranded nucleic acid target molecule havinga sequence and a plurality of binding sites disposed along the sequence.A plurality of probe molecules having a first sequence specificity areadded to the double-stranded nucleic acid target molecule. The probemolecules having the first sequence specificity and the target moleculeare incubated so as to effectuate preferential binding of the firstprobe molecules to both a first binding site and a second binding siteof the target molecule. A parameter related to a distance between thefirst binding site and the second binding site is measured.

A second probe molecule may be added to the double-stranded nucleic acidtarget molecule, and the second probe molecule and the target moleculemay be incubated so as to effectuate preferential binding of the secondprobe molecule to both a third binding site and a fourth binding site ofthe target molecule. A parameter related to the distance between thethird binding site and the fourth binding site may be measured. Theparameter may be measured by passing the target molecule and bound probemolecules through a nanopore.

The description of elements of the embodiments of other aspects of theinvention may be applied to this aspect of the invention as well.

BRIEF DESCRIPTION OF THE DRAWINGS

The foregoing features of the invention will be more readily understoodby reference to the following detailed description, taken with referenceto the accompanying drawings, in which:

FIG. 1 is a flow chart for a method of sequencing DNA in accordance withan embodiment of the present invention;

FIG. 2 is a schematic drawing showing a double-stranded DNA targetcomplexed with multiple probe molecules to form local ternary complexes;

FIG. 3 is a schematic drawing showing a double-stranded DNA targetcomplexed with multiple strand-invading probe molecules to form localternary complexes;

FIG. 4 a is a schematic drawing showing an idealized kinetic trace of animpedance measurement generated as local probe-target complexes passesthrough a nanopore;

FIGS. 4 b and 4 c each show a graph of an exemplary current signaldetermined in accordance with an embodiment of the invention;

FIG. 5 a is a schematic diagram of a nanopore detector in accordancewith an embodiment of the invention;

FIG. 5 b is a schematic diagram of the detector of FIG. 5 a with a localcomplex within the nanopore in accordance with an embodiment of theinvention;

FIG. 5 c is a schematic diagram of the detector of FIG. 5 a with anuncomplexed region of target DNA within the nanopore in accordance withan embodiment of the invention;

FIG. 6 is a flow chart illustrating a reconstruction algorithm inaccordance with an embodiment of the present invention; and

FIG. 7 is a flow chart illustrating an algorithm for reducingpossibilities during DNA sequencing in accordance with an embodiment ofthe present invention.

DETAILED DESCRIPTION

As used in this description and the accompanying claims, the followingterms shall have the meanings indicated, unless the context otherwiserequires:

A “target” means a biopolymer, for example, having sequence informationwhich is to be determined using embodiments of the present invention;

A “probe” means any molecule or assembly of molecules capable ofsequence-specific covalent or non-covalent binding to a target molecule;

In connection with a target bound to one or more probes, a “ternarycomplex” means a complex formed from three or more complementarybiomolecules. For example, a complex of two complementary strands of DNA(i.e., double-stranded DNA or dsDNA) with a third probe strand is aternary complex. Higher order complexes also include ternary complexes;for example, a complex of double-stranded target DNA with two probemolecules bound to the same stretch of sequence, with two probemolecules bound to each other and in turn to a stretch of sequence, orwith two identical probe molecules bound to different targetsubsequences are quaternary complexes that also include a ternarycomplex.

A “local ternary complex” means a ternary complex that is localized to aparticular stretch of biopolymer. For example, multiple local ternarycomplexes may exist when probe molecules bind to a target molecule inmultiple locations separated by stretches of bases to which the probe isnot bound (i.e., uncomplexed regions).

A “polyamide” refers to the class of DNA binding polyamides originallydeveloped by the Dervan laboratory at the California Institute ofTechnology. Polyamides are molecules containing heterocycle ringstructures that can be combined in a modular fashion to recognize DNAsequences by binding in the minor groove of duplex DNA. Typicalheterocycle ring structures in this class are, but not limited to,N-methylimidazole, N-methylpyrrole, 3-hydroxy-N-methylpyrrole,N-methylpyrazole, 3-methylthiophene, benzimidazole,hydroxybenzimidazole, imidazopyridine. The rings may be connected viacarboxamide linkages, by the amino acid β-alanine, or directly joined bysingle bonds. Other ring structures and linkages can be envisaged bypractitioners familiar with the art.

The polyamides may include spacers that are not sequence specific andthat allow the curvature of the polyamide to stay in register with thecurvature of double-stranded DNA (“gapped polyamides”). The spacer maybe used as a gap when the polyamide is used as a probe for sequencedetermination. Gapped probes are more efficient at reconstructingsequence information than a non-gapped probe.

In connection with a target and a probe, a “probe map” means a data setcontaining information related to the sites along a target sequence atwhich a probe preferentially binds. The data set may include absolutepositional information referenced to a known sequence, relativeinformation related to distances between binding sites, or both. Thedata set may be stored in computer media.

A “spectrum map” means a collection of probe maps determined for atarget sequence.

A probe recognition site refers to a target sequence, structure, orconformation that is preferentially recognized by the probe in thepresence of other target sequences, structures, or conformations. Theprobe may interact with the recognition site through non-covalent,covalent, or mixtures of covalent and non-covalent interactions.

In connection with a probe map or a spectrum map, “dynamically weightedspectrum” means the collection of probes under consideration for theextension of a growing sequence with weights as to the likelihood thatthey should be chosen next. These weights will be based on a parameterrelated to the measured positions of the probes. The dynamicallyweighted spectrum may or may not include sequence information. Forexample, in the case of six-mers, if the growing sequence ends in ATACG,the dynamically weighted spectrum may include only the probes ATACGA,ATACGC, ATACGG, and ATACGT with weights based on their (relative)positions. Alternatively, the dynamically weighted spectrum may bethought of as a collection of probes, independent of sequence, whoseprobability of being next in the sequence is based on their measureddistances from the current growing end of the strand. (In this casesequence-specific information would be incorporated subsequently.)

Illustrative embodiments described herein relate to detecting andanalyzing probe maps, to gather sequence information. More particularly,illustrative embodiments described herein relate to creating anddetecting specific ternary complexes along double-stranded biopolymertarget molecules to gather sequence information. For simplicity, theembodiments described below use DNA as a target, but other biopolymerssuch as RNA may also be sequenced. In an embodiment, probes that areselective for one or more sequences are hybridized to a target sequence.The target sequence may be double-stranded DNA and the probes are chosenfor the ability to bind to one or more sites of double-stranded DNA. Thebinding of the probe is then detected using a nanopore to generate aprobe map in computer media. By detecting probe binding at multiplesubsequences along the target, a spectrum map may be constructed.Information contained in the spectrum map may then be aligned anassembled into a finished sequence.

Like SBH, positional hybridization sequencing relies on the biochemicalhybridization of probes to an unknown target. However, once thebiochemical hybridization is substantially complete, a detector thatdetermines the position of hybridization is utilized. In addition todetermining the position of hybridization, the detector is capable ofcounting the number of times a probe hybridizes to the target. Thepositional information allows a combinatorial sequence reconstructionalgorithm to be used that differs from those used in SBH. In standardSBH, the spectrum of the target sequence from which probes are selectedduring each step of sequence reconstruction includes all the probes thathybridized and this remains constant during the course of thereconstruction. That is, every probe in the spectrum can be used at anystep during reconstruction. In the case of positional hybridizationdetection, the spectrum is a dynamically weighted spectrum. Thedynamically weighted spectrum may include all the probes that bound tothe target (with different weights). However, in most embodiments thedynamically weighted spectrum will only include a subset of thespectrum. A consequence of not including all the probes is that duringreconstruction of the sequence the dynamically weighted spectrum mustchange at least once in order to encompass all probes in the spectrum.

Other illustrative embodiments relate to methods and computer algorithmsfor aggregating, aligning and assembling multiple probe maps into amodel nucleotide sequence. In a specific embodiment, combinatorialsequence reconstruction algorithm is used to generate spectrum maps fordetermining sequence information. Due to the complex nature of thecalculations and quantity of data generated, it is expected that themethods described herein be carried out in a substantially automatedmanner by an appropriate computer algorithm running on a computer ofadequate speed.

In an illustrative embodiment, detection of local ternary complexes isaccomplished by passing double-stranded biopolymer target molecule orfragment thereof through a nanopore and detecting an electrical signalindicative of the locations of the local ternary complexes along thedouble-stranded biopolymer target molecule. The use of nanopores todetect solution phase DNA hybridization reactions is described in U.S.Pat. No. 6,537,755 to Drmanac, and in U.S. Patent ApplicationPublication No. 20060287833, to Yakhini, both hereby incorporated hereinby reference in their entirety.

FIG. 1 shows a flow chart for a method of sequencing a target inaccordance with an embodiment of the present invention. For simplicity,this and other embodiments detailed herein use genomic DNA as the targetto be sequenced, although other types of molecules, including cDNA, andRNA could be sequenced in this way. A target DNA sequence is extractedfrom a source and purified (step 100). The starting material may be anyof a variety of single or double-stranded nucleic acids, for example,genomic DNA, PCR products, cDNA, RNA-DNA hybrids and the like. The useof double-stranded DNA as a target has several advantages over the useof single stranded DNA. These include a reduction in the amount ofsecondary structure present in the target strand to be sequenced, a moreconsistent transmission speed of the DNA through the nanopore, theability to manipulate single molecules of double-stranded target DNAprior to introduction into the nanopore, a longer persistence length ofthe DNA, and the use of probe molecules that are of higher affinity andsequence selectivity. While the method may be most powerful with longpieces of DNA to which a probe may bind in multiple subsequencepositions, the invention is also applicable to shorter nucleic acidssuch as viral genomes, plasmids, cosmids etc. The target DNA may bederived from any of multiple sources and may be the genome of anorganism whose DNA sequence is unknown, or is known (resequencing).Alternatively, the DNA may be derived from several organisms. The targetDNA may be of a known size or range of sizes; alternatively, the size ofthe target DNA may be unknown. The purification of step 100 may be thecomplete purification step or performed in conjunction with later steps.For example, a biochip that includes a nanopore may also accomplishfurther purification via electrophoretic or other methods. Optionally,the target may be amplified by PCR, rolling circle amplification, orother suitable method.

If chromosomal DNA or other large target DNA polymers are to besequenced, it may be desirable to fragment the DNA (step 120). TargetDNA may be fragmented by any of a number of commonly used methodsincluding hydrodynamic shearing, sonication, ultrasonic fragmentation,enzymatic cleavage, nebulization, chemical cleavage, and heat inducedfragmentation. Hydrodynamic shearing may be favorable if large nucleicacid fragments are desired. A commercially available device, HydroShear(Genomic Solutions, Ann Arbor, Mich., USA), is available that will shearDNA to a tight size distribution. To use this device, double-strandedDNA in solution is passed through a tube with an abrupt constriction.Fluid accelerates to maintain the volumetric flow rate through thesmaller area of the constriction. During this acceleration, drag forcesstretch the DNA until it snaps. The DNA fragments until the pieces aretoo short for the shearing forces to break the chemical bonds. The flowrate of the fluid and the size of the constriction determine the finalDNA fragment sizes. In one configuration of the device, size ranges ofsheared DNA of from 650 base pairs to 40,000 base pairs are achievable.In a second configuration of the device, size ranges of 1,000 base pairsto 9,000 base pairs are achievable.

The target DNA may then be aliquoted (step 120) and combined with probemolecules (step 130). For example, a first aliquot may be combined witha solution containing molecules of probe specific for a first targetsubsequence (e.g., a 5 to 20 base pair subsequence) and a second aliquotmay be combined with a probe that is specific for a second targetsubsequence. In a simple scheme, the number of aliquots is equal to thenumber of different probes to be used in the experiment. A single probe(i.e., a collection of probe molecules having the same or functionallyequivalent sequence specificity) may then be added to each aliquot ofthe target DNA. Alternatively, the number of aliquots may be smaller orlarger than the number of different probes; a mixture of probes may beadded to each aliquot. The probes may be of either uniform or differinglengths. Various mixing procedures may be used to combine the target andprobe, but if the procedure is performed at very low volumes (e.g.,0.1-10 nanoliter), diffusional mixing may suffice.

The aliquots are then incubated under stringent conditions (step 140) toallow the probe molecules to sample multiple subsequences and toapproach equilibrium in finding energetically favored subsequences tobind to (i.e., the subsequences for which the probe is specific). Inother words, each aliquot containing target DNA and a probe or a mixtureof probes is incubated under conditions that result in greaterhybridization of probes to matched recognition sequences than tosequences containing a mismatch. Taking into account the type of probeused, stringency may be adjusted through temperature, saltconcentrations, addition of organic solvents, washing with solutions,electrophoretic washing, or other methods. Optionally, incubatedaliquots may be combined prior to analysis, especially if tagged probesare used. The result is a formation of specific local probe-targetcomplexes, which, as described below, may be modeled as a probe map. Inspecific embodiments, the probe map comprises the relative positionalinformation of the local ternary complexes along the double-strandedbiopolymer, the absolute positional information of the local ternarycomplexes along the double-stranded biopolymer, and error of positionalinformation of the local ternary complexes along the double-strandedbiopolymer. In another embodiment, a second probe specific for one ormore second recognition sites of the target molecule is used to generatea second probe map. In yet another embodiment, two or more probes areused for the detection of local ternary complexes, and for thegeneration of two or more probe maps.

The collection of hybridized or otherwise bound probes and positions maybe used to form a model probe map of the target DNA. The probe map maybe used to reconstruct the target DNA sequence. For example, multipleprobe maps from probes having a variety of sequence specificities may bealigned and assembled to reconstruct a sequence. In an embodiment, n-merprobes are used (e.g. polynucleotide probes), each probe having aspecificity for one of 4^(n) possible DNA sequence combinations. Acomplete library of 4^(n) n-mer probes may be utilized to give completetarget sequence coverage. In illustrative embodiments, n is from 3 to10.

Since the target was fragmented, the complete probe map will beconstructed from data derived from the multiple fragments. The probe mapmay be assembled by matching the pattern between overlapping regions ofmultiple fragments to arrive at an aggregated decoration pattern for agiven probe having a given target specificity. The pattern matchingroutine may include determining the 3′→5′ orientation of each decorationpattern, especially if this has not been determined previously, e.g. bycomparison to a known sequence in a resequencing operation.

The sequence specificity of a probe, combined with its estimated probemap with respect to the target yields a partial sequence model of thetarget. In other words, if a probe-target complex is detected, it may beexpected that the target sequence contains the subsequence for which agiven probe is specific. For sufficiently short probes, this subsequencewill likely be repeated multiple times in the target at positionscorresponding to the positional information that should be availablefrom detection of binding of that probe. The subsequences may beassigned relative positional information by making use of the distancesof measured uncomplexed regions intervening between the probe bindingsites. The positional information will likely contain some errorhowever.

The model sequences derived from the probe map of multiple probes havingmultiple corresponding sequence specificities may be preliminarilyaligned with respect to each other and/or with respect to otherlandmarks (step 150). For example, for resequencing operations, thesequences may be aligned based on a maximal overlap with a knownsequence. Alternatively, the model sequences may be aligned based onmaximal overlap with each other or with respect to bound proteins (e.g.zinc finger proteins) included in the mixture, or with respect to othersequence landmarks.

A starting model sequence is created (step 160). For example, a startingsequence may be created using the known sequence specificity of one ofthe probes for which a probe map has been determined. The chosen probemay be located at or near the 3′ or 5′ terminus of the preliminarilyaligned spectrum map. A growing-end may be chosen for the startingsequence, on the internal side.

In certain embodiments, a spectrum map is generated from at least afirst probe map and a second probe map. Determining the spectrum map mayinclude using a combinatorial sequence reconstruction algorithm.Determining the spectrum map may also include determining a candidatesequence by ordering at least two probe sequences using positionalinformation and error of positional information. If the firstrecognition site sequence and the second recognition site sequenceoverlap, determining the spectrum map may include aligning therecognition site sequences.

A weighted spectrum is chosen as a subset of the spectrum map withrespect to the growing end (step 170). The probes encompassed in thissubset of the spectrum are selected using positional information or someparameter related to positional information generated during thedetection step. The various sequences corresponding to theprobe-specificities for the probes for which binding events weredetected within the weighted spectrum are then compared with a sequenceof the growing end (in the initial stages of the process, this may bethe entire starting sequence). If possible, a sequence is chosen thathas maximal identity for the growing end yet extends at least onenucleotide beyond the growing end (step 180). At least one base of thenon-overlapping sequence is then appended to the growing end (step 190).The weighted spectrum may then be changed. This process of selectingnucleotides, appending nucleotides, and changing the weighted spectrumis repeated until sufficient sequence is determined or the data set isexhausted. If more than one nucleotide choice is available forappending, this may be resolved by evaluating each possible choice interms of maximal overlap for the next few sequence extensions.

Contrastingly, in standard SBH, the weighted spectrum of the targetsequence from which probes are selected during each step of sequencereconstruction is exactly equal to the spectrum of the target. That is,every probe in the spectrum can be used at any step duringreconstruction. In the case of positional hybridization detection, theweighted spectrum of the target sequence encompassing probes selectedduring each step of reconstruction may be smaller than the length of thetarget strand. The weighted spectrum may be as large as the spectrum ofthe target minus one probe. However, in most embodiments the weightedspectrum will be significantly smaller than the target spectrum. Aconsequence is that during reconstruction of the sequence, the weightedspectrum must change at least once in order to encompass all probes inthe spectrum.

If there were no error in the probe maps with regard to the position ofhybridization of the probe, the weighted spectrum could consist of onlyone probe. Some error is to be expected, however. At the other extreme,if there were no certainty regarding the position of probe locations ina probe map, then the next nucleotide to be appended would logically beselected from any probe binding to the target. This would correspond tostandard SBH. Balancing these two extremes, an embodiment of theinvention optimizes the weighted spectrum based on estimated error inthe positional information contained in each probe map or in thespectrum map. For example, a standard deviation, variance or otherquantitative measurement of the error in one or more probe locations maybe calculated based on the detector output. The weighted spectrum maythen be set to be a multiple of this estimated error. Accordingly, thenumber of probes encompassed by the weighted spectrum and the size ofthe error may be inversely proportional, or otherwise inverselycorrelated. Simply stated, a larger weighted spectrum may be used when alarger positional error is present. The weighted spectrum may beexpressed in terms of probes that fall within a fixed number ofnucleotides distance from the growing end. Alternatively, the range mayvary based on the error in the positional information for each probe, oreven based on the estimated error at each probe binding site.

Other alternative calculation methods may be employed including adaptiveor genetic algorithms that adjust based on the quality of the sequenceoutput. The output of multiple such methods may also be compared and thebest one selected. The range or ranges may be calculated once, multipletimes during the procedure, or after each nucleotide is appended.

A number of types of molecules are known to bind in a sequence-specificfashion to double-stranded DNA and may be used as probes in the currentinvention individually or in combination. The molecules may be composedof one or more moieties such as polyamides, a nucleic acid, and/or apeptide nucleic acid. In a specific embodiment, the molecule is a gappedpolyamide molecule such as those described in the following Dervanreferences: Dervan, P. B. Molecular Recognition of DNA by SmallMolecules. Bioorg. Med. Chem. 2001, 9, 2215-2235; and Dervan, P. B.;Edelson, B. S. Recognition of the DNA minor groove by pyrrole-imidazolepolyamides. Curr. Op. Struc. Biol. 2003, 13, 284-299.

Alternatively, the probe molecules may be composed of native or modifiedDNA oligonucleotides that are capable of binding in a triple strandedfashion with double-stranded DNA. Alternatively, the probe molecules maybe nucleic acid molecules that are combined with the target in thepresence of a recombinational enzyme or a recombinational protein. Anon-limiting example of a useful recombinational enzyme is recA. Furtherinformation about binding recA probes to DNA strands is provided in thefollowing two references, incorporated by reference in their entiretyherein: Bianco et al., Nucleic Acids Research Vol. 24, No. 24 (1996)4933-4939; and Cox, Nature Reviews—Molecular Cell Biology, Vol. 9(February 2007) 127-138. Alternatively, the molecule may be a peptidenucleic acid capable of binding to double-stranded DNA. Other modifiedoligonucleotides, proteins, peptides, or other polymers are also capableof binding to double-stranded DNA in a sequence specific fashion andthus may be used.

As shown schematically in FIG. 2, probes 200 will bind to specificsubsequences along a DNA strand 210 to create a probe map. In anembodiment, the DNA is double-stranded DNA and the probes may bind toboth strands with comparable affinity or preferably to one strand, andmay bind to the minor or major groove of helical dsDNA targets.Additionally, as shown in FIG. 3, the probes may invade and disrupt thebase pairing of the dsDNA, as is known to occur with certain peptidenucleic acids. In an embodiment, the probes may contain moieties that donot have sequence specific binding behavior. The probe may contain twosequence-specific moieties separated by a relatively non-specific linkermoiety. In the case of an oligonucleotide, the non-specific bindingregion may be composed of a universal base that is able to bind to allfour cognate bases (see, e.g., Preparata U.S. Patent Application No.20010004728.) In the case that the target is double-stranded DNA, theprobes may bind to both strands with comparable affinity or preferablyto one strand, and may bind to the minor or major groove of helicaldsDNA targets. Additionally, as shown in FIG. 3, the probes may invadeand disrupt the base pairing of the dsDNA, as is known to occur withcertain peptide nucleic acids. Non-specific linkers may also be includedin probes that include binding regions composed of polyamides, peptidenucleic acids, peptides, or proteins. Probes may also be labeled toincrease detection sensitivity, or to allow discrimination of differentprobes having different corresponding sequence specificities.

In an embodiment, the probes are polyamides that bind in a 2:1 complexwith dsDNA in the minor groove of the DNA. Alternatively or in addition,the polyamide may bind to the DNA in a 1:1 complex or in a 2:1 orhigher-order complex. In the case of a 2:1 complex, the two polyamidesmay be separate molecules or may be covalently joined at the ends or inthe middle of the molecules to form hairpin polyamides, cyclicpolyamides, H-pin motifs, or U-pin motifs. Charged groups may be placedat either end or in the middle of the polyamides in order to changetheir affinity for target or for other probe molecules. Linkers may beattached between heterocycles in the polyamide to modulate the curvatureof the molecule. Curvature modulation may be used to increase thebinding affinity of a polyamide to a longer DNA sequence by matching thecurvature of the polyamide to the curvature of the DNA molecule. Otherlinkages and appendages may also be used.

Polyamide probes may consist of imidazole, pyrrole, hydroxy-pyrrole,benzimidazole, and hydroxybenzimidazole heterocycles. The heterocyclesmay be substituted with alkyl or functional groups in order to modulatetheir binding to DNA. Other heterocycles may also be used. Pairing rulesof polyamides for minor groove recognition of DNA have been described byDervan. (Dervan, P. B.; Edelson, B. S. Recognition of the DNA minorgroove by pyrrole-imidazole polyamides. Curr. Op. Struc. Biol. 2003, 13,284-299.)

Generally, each aliquot of target DNA is subsequently analyzed to detectbinding of the probe to the target DNA. In a specific embodiment, aprobe map is generated for each probe to give absolute or relativepositional information of probe binding sites located on a longerstretch of dsDNA. The probe binding sites are separated by non-binding,uncomplexed regions of the target that have relatively low affinity forthe probe. Accordingly, the probe map may include absolute or relativeinformation related to the location and length (e.g., in base pairs) ofthe uncomplexed regions.

In an embodiment of the invention, a single molecule detector is used todetect probe binding and thus arrive at a probe map. Unlike detectorstypically used for SBH, the single molecule detector may determine thenumber of times that a given probe binds to a target molecule and theabsolute or relative binding sites for that probe. For example, ananopore based detector, a field effect transistor, Coulomb-chargingbased detector, or scanning-tunneling or force microscope may be used.Multiple such detectors may be used in parallel or in series. Data frommultiple detectors may then be aggregated for analysis, includingcorrelation analysis. Each mixture may be independently analyzed. In anembodiment of the technique, a nanopore is used to detect the positionof local ternary complexes along on the target DNA.

In certain embodiments, a nanopore-based detector detects temporalchanges in current as probe-complexed and uncomplexed stretches ofdouble-stranded DNA pass through the pore. As seen for single strandedDNA, the change in current may be positive or negative depending on theconcentration of ions in the buffer on either side of the nanopore. Thecurrent may increase or decrease while a probe is passing through thepore due to blockade current, tunneling current or other mechanism.

FIG. 4 a shows an idealized trace of a nanopore readout in terms ofimpedance as a function of time. Each time a probe molecule 200 passesthrough the nanopore detector, the impedance level increases (i.e.,current blockade). Accordingly, the nanopore signal may be elevated ordepressed for a period of time that reflects the length of theprobe-target complex (item 400), or the length of the interveningnon-binding regions (item 410). Relevant readout information may berecorded on computer media for future analysis.

A current signal, reflecting a change in impedance or current measuredby the nanopore detector, may be monitored to gain information about thepresence and spacings of hybridized probes. As an analyte moves througha volume monitored by the nanopore detector, the current signal changes.The signal may be elevated or depressed for a period of time thatreflects the length of the analyte, e.g., a probe-target complex, or thelength of the intervening regions without probes. A typical analyte willimpede the flow of ions in the electrolyte and is non-conductive.Therefore, the current typically decreases as the analyte flows throughthe sensing volume. In some embodiments, e.g., a low salt electrolyteand a charge-carrying analyte, the current signal may increase as theanalyte flows through the sensing volume. The current signal furtherchanges when the portion of the analyte containing the hybridized probemoves through the volume between the sensing electrodes.

Referring to FIG. 4 b, sample data is presented for a piece ofdouble-stranded DNA with an attached probe passing through a pore. Thex-axis is time in milliseconds and the vertical axis is current inpicoamps. A negative current was run, so that the trace starts at about−3920 pA. That is the current when only buffer is present in the port.As a DNA strand entered the pore, the current was reduced by almost 200pA to about −3720 pA. When the section of DNA that had the probeattached passed through the nanopore, there was a further reduction inthe current to a little less than −3420 pA. After the DNA and probepassed through, only DNA was in the nanopore, and the current returnedto −3720 pA. When the DNA exited the pore, the current returned to thebaseline of −3920 pA. FIG. 4 b provides similar data, with theapplication of a baseline correction. Thus, the baseline appears at 0 pAand the peaks indicate the difference between the baseline and thechosen peak height.

A time interval between current signal changes may be recorded. Theduration of the change in the current signal may indicate a presence ofa hybridized probe. This duration may be used to determine a distancebetween two probes on the biopolymer. To determine the distance, one maycalibrate the system with known standards and calculate an average speedfor the analyte in the channel. Time would then be directly related todistance. In embodiments with multiple hybridized probes, multiple peaksmay form, with each peak corresponding to a hybridized probe, with thetime between peaks being indicative of the relative spacing of thehybridized probes.

Similarly, the duration of a change in the current signal may be used todetermine a length of the analyte.

FIGS. 5 a-5 c schematically show a nanopore detector in accordance withan embodiment of the invention. A cathode 500 and anode 510 (e.g.,platinum terminals connected to an appropriate power supply) arepositioned to create an electrophoretic field in a buffer solution. Thesolution is divided into two chambers by a nanopore 530. As a decoratedtarget-probe complex 520 is electrophoretically driven through thenanopore 530 by the electrophoretic field, a detection circuit 540detects and records positional information related to the location ofthe bound probes along the length of the target to give a probespectrum. The target-probe complex is shown as double-stranded, butcould be the single stranded complex of FIG. 2. FIG. 5 a shows thecomplex 520 prior to entering the pore 530. FIG. 5 b shows the complex520 within the pore 530; in this state the detection circuit 540 shoulddetect the presence of a bound probe (through an increase or decrease ofsignal). FIG. 5 c shows the complex in a configuration for which thedetection circuit 540 should detect the absence of a bound probe. Theorientation of the target as it passes through a nanopore may bedetermined from the probe map, e.g. by comparing the map to a knownsequence in a resequencing procedure.

In a related embodiment, the target-probe complex 520 passes throughmultiple “stacked” nanopore detectors in sequence. In this way, multipledata sets may be collected and analyzed to give a better signal to noiseratio. This may result in greater accuracy in the probe map.Alternatively, other single molecule detection systems may provide suchcorrelative data.

In a certain embodiment, the biopolymer target molecule is singlestranded DNA (ssDNA). Double-stranded DNA may be generated from thesingle-stranded template by enzymatically synthesizing the second strandand thereby afford may of the aforementioned benefits of using dsDNA: areduction in the amount of secondary structure present in the target, amore consistent transmission speed of the DNA through the nanopore, theability to manipulate single molecules of target DNA prior tointroduction into the nanopore, and a longer persistence length of theDNA. In an alternative embodiment, ssDNA is stabilized using a pluralityof single-stranded DNA binding agents to stabilize the single strandedDNA. The binding agent molecules may be nonspecific, or may be chosen tohave specificity for various DNA target subsequences. For example, amixture of tens, hundreds, or thousands of short (e.g., designed to bindto 6-mers) DNA, or DNA analogs designed to target particular targetsubsequences may be stringently incubated with a target ssDNA sample.One or more detection probes may then be added and the mixture incubatedunder stringent conditions. The mixture may be chosen so as to notpreclude a detection probe from binding to its complementarysubsequence. Different stabilizing agent mixtures may be used withadditional aliquots containing target combined with additional detectionprobes having different subsequence specificity. In related embodiments,the mixtures of stabilizing agents may be chosen to be non-overlappingin their sequence specificity. Using locked nucleic acid or peptidenucleic acid probes may allow the use of longer and more stable probes,capable of binding in higher stringency conditions (e.g., hightemperature or low-salt). As a result, the mixture may require feweragents; e.g., a mixture may be sufficient that binds to only a third orless of the ssDNA sequence.

The reconstruction algorithm, in accordance with some embodiments of theinvention, may include the following steps. The reconstruction algorithmproceeds iteratively, adding letters to the reconstructed DNA sequencefrom 5′ to 3′. Progress along the sequence is measured by the variableposition, which increases in increments of 1 from 0. The set ofavailable probes is the population of probes whose annealing to thetarget strand is measured at a position within some distance Δ ofposition. Each element in the set of available probes is a Record, anobject containing the probe's sequence and the position at which it wasmeasured on the target strand. Let Possibility denote a structureconsisting of a string denoting a DNA sequence, a real number scoreproportional to the logarithm of Possibility's probability, a booleanvalue sentinel, and a list of Records used by this sequence.

Referring to FIG. 6, the algorithm 600 proceeds in three stages:initiation, elongation and termination.

Initiation:

Position is initialized to 0 and a list of Possibilities with onePossibility per Record in the set of available probes about 0 (step610). The Possibility with Record r is given score of (r.position−0)^2.

Elongation:

Increment position by 1 (step 620).

Trim the list of Possibilities as described below with reference to FIG.7 (step 630).

For each Possibility p in the list of Possibilities (step 640):

-   -   Let S denote the set of Records r in the set of available probes        that satisfy the condition that, if l is the length of the probe        indicated by r, the l−1-length suffix of the sequence denoted by        p is equal to the l−1-length prefix of the probe indicated by r        (step 650). (Note that it is not stipulated that all probes are        of a length, as discussed in more detail below.) Such a probe        sequence is called an extension of p. r is also preferably        absent from p's list of Records used, so that each Possibility        uses each Record only once. If more than 1 Record in S denotes        the same probe-sequence, remove all but the one with the lowest        measured position.    -   If |S|=1 (step 660), let the single Record in S be called r.        Append the final character of r's probe to p's DNA sequence, add        (r.position−position)^2 to p.score, set p.sentinel to false and        add r to the set of Records used by p (step 665).    -   If |S|=0 (step 670) and p.sentinel=true, remove p from the list        of Possibilities and continue. If p.sentinel=false, remove p        from the list of Possibilities and add 4 new Possibilities to        the list with the DNA sequence and score ofp, setting their        sentinels to true. Increase the score of each new Possibility by        some fixed amount (Δ^2 may be used). This action (extending a        sequence without finding a Record to extend it) is called a        reach (step 675).    -   If |S|>1 (step 680), remove p from the list of Possibilities.        For each Record r in S, add a new Possibility q which is a copy        of p. Append the final character of r's probe to q's DNA        sequence, add (r.position−position)^2 to q.score, set q.sentinel        to true and add r to the set of Records used by q. This action        (extending new sequences according to multiple found Records to        extend it) is called a branch (step 685).        Termination:

Eventually, the end of the data set is reached and the set of availableprobes will be empty (step 690). At this point, the list ofPossibilities will decrease in membership until there are noPossibilities left. As Possibilities are removed, a reference is kept tothe best Possibility. best is the Possibility whose number of usedprobes is greatest, measured by position at which the Possibility isterminated minus the number of reaches the Possibility made. If morethan one Possibility have used the same number of probes, best is thePossibility with the lowest score. The sequence stored in best isreturned as the reconstructed sequence of the target strand (step 695).

Reducing List Size

As described above, the algorithm 600 has exponential runtime becausethe list of Possibilities grows by a factor of at least one at eachposition. In practice, the base of this exponent is highly sequence andprobe-length-dependent since a branch takes place when two probes withthe same prefix are found within a distance of 2Δ. Since exponentialruntimes are intractable, it is preferable to keep the size of the listof Possibilities from growing too much. Referring to FIG. 7, at eachstep, before elongating the putative Possibilities (step 700), thefollowing list-size reducing measures (step 630) are enacted:

A. Haircutting (step 710)

-   -   Since a Possibility's score is proportional to the        log-probability of the Possibility, one may impose a lowest        probability that one is willing to entertain in a Possibility.        All Possibilities in the list whose probabilities are worse than        the best-scoring Possibility by a certain amount are removed        from the list (step 720).

B. Parallel culling (step 730)

-   -   A problem with traditional SBH using probes of length k are the        so-called SBH failure modes. These modes arise in sequences of        the form AbAcA or AcBdAeB, where A and B are sequences of length        k−1 and b, c, d and e are sequences of any length. In the first        case, it is not possible to distinguish between AbAcA and AcAbA        while in the second case, it is not possible to distinguish        between AcBdAeB and AeBdAcB. These modes pose a challenge to the        disclosed algorithm as well, because each SBH failure mode        produces two Possibilities that have used exactly the same        probes and now extend the same sequence. Therefore, each SBH        failure mode doubles the number of Possibilities in the list. In        parallel culling, the complete list of active Possibilities is        examined. Any two Possibilities that have used exactly the same        Records and are extending the same sequence (i.e. S will be        identical when elongation is applied to these Possibilities),        are a case of parallel paths. The future, including all reaches        and branches of these paths are the same, with the exception        that their scores differ by a constant amount. As such, one has        all the information one is going to get about the difference in        probability between these paths, and one preferably picks the        better-scoring option (step 735) and discards the other (step        720). The code is also outfitted to search a reference human        genome and eliminate the Possibility whose most recent sequence        extensions are not found in the reference. This is effective        because, while existing sequencing technologies make large-scale        errors, they are generally accurate on the scale in which SBH        failure-modes arise in algorithm 600.

C. Jump-forward culling (step 740)

-   -   Consider two Possibilities p and q that diverged at a branch.        Without loss of generality, assume that p is the correct path        and that q has been extending incorrect sequence since this        branch. Since the re-use of Records is not permitted, the first        Record used by q and not by p preferably comes from later on in        the sequence and not earlier on. As such, all branches away from        the correct path are jump-ahead errors. Let the position at the        branch be t and the true position of the first Record used by q        after the branch be t+u (true position denotes the position at        which the correct path would incorporate the Record). The u        probes in between t and t+u remain untouched by q unless q        enacts a jump backward to use them. Note that if q makes this        jump back and a subsequent jump forward to use all the probes        used by p, q and p will have used the same set of probes and one        will be in an SBH-failure-mode situation, already handled above.        Thus, assume that q does not jump back for these u Records. As        such, at any position, q is using a probe whose true position is        greater than that of the probe used by p.        -   Once position increases to the point where all u of the            Records skipped by q have positions too low to be in the set            of available probes, the set of probes used by q that are in            the set of available probes will be a superset of the probes            used by p in the set of available probes. Another way of            thinking about this is that because q has skipped over u            Records and continues to extend the same number of probes            asp, q has u fewer Records in its future. The set of Records            in the future is the union of the Records whose positions            are too high to be in the set of available probes and the            set of unused Records in the set of available probes. So            having a greater number of used Records inside the set of            available probes is identical to having fewer Records            available for future use. Since both p and q have used the            same number of probes (they are both extending sequences            position long) and p has more Records available for future            use, p will use more Records in the total and extends a more            correct sequence than does q. Therefore, it is safe to            remove all Possibilities q from the list for which there            exist a Possibility p in the list that uses a subset of the            Records still in the set of available possibilities that q            does.            Variable Probe-Lengths

In traditional SBH, all probes are of a common length, usually called k.Here, this constraint has been removed as a result of observations aboutthe nature of human DNA. Simulations of DNA sequencing algorithms onrandom DNA (data strings in which each letter has an independentprobability of being “a,” “c,” “g,” or “t”) are inevitably morefavorable than simulations done on sequences drawn from human DNA. Inthe case of the disclosed algorithm, this is because short sequences arefound repeated in very short stretches. In classical SBH, a repeatedelement of length k−1 produces a branch-point. In the disclosedapproach, this increases the computational burden by doubling the numberof Possibilities in the list. It also exposes one to the chance ofmaking incorrect decisions in cutting down the size of the list byHaircutting or Parallel-culling. In random DNA, the probability offinding two Records in the set of available probes with the samek−1-length prefix is equal for all k−1-length prefix strings. In humanDNA, however, certain strings are enormously more common. To avoid theseambiguities, longer probes are used for all sequences exhibiting higherfrequency.

A. Non-repeated prefixes

-   -   Consider the case where the input string has the following form:        acatagtctgat (SEQ ID NO: 1). gtatagtctgaa (SEQ ID NO: 2) (where        . . . is a sequence of DNA no longer than A). One may wish to        sequence this string by classical SBH with k=8. The underlined        prefix agtctga is found twice, which would introduce a branch in        the reconstruction. To avoid this ambiguity, one may extend the        probe prefix back until the two are different: (agtctga,        agtctga)->(catagtctga (SEQ ID NO: 3), tatagtctga (SEQ ID NO:        4)). Using a probe of length 11 instead of 8 eliminates this        branch point and restores unambiguous reconstruction of this        sequence. By analyzing existing versions of the human genome,        one can develop a collection of proximally-repeated short        elements (“ccctccct” is an identified example, as are all of the        homo-polymer strings) for which to recommend longer probes.

B. Information Content

-   -   Another approach for probe set design is an examination of the        information content of certain probes. In the human genome,        there are functionally-relevant regions of reduced alphabet,        i.e., regions in which only 2 or 3 of the DNA bases are found.        These regions are generally either structural DNA (centromeres,        telomeres, constitutive heterochromatin) or DNA that codes for        functional RNA. In such regions, all probes exhibit the reduced        alphabet of the region, which makes for greater ambiguity. For        example, there are 4^7=16,384 prefixes of length 7 in a region        in which the alphabet is all four letters of the DNA alphabet        but only 2^7=128 prefixes of length 7 for a region with alphabet        reduced to size 2. Therefore in a region of alphabet of size 2,        there are only 128 possible prefixes and a branch event is 128        times as likely as in a 4-alphabet region. To increase the        number of prefixes, it is proposed that probes with reduced        alphabet be longer. Specifically, if full-alphabet probes of        length 8 were used (prefix of length 7), 3-letter-alphabet        probes of length 10 would be used (3^9=19,683˜16,384) and        2-letter-alphabet probes of length 15 would be used        (2^14=16,384).

C. Non-Redundant Prefixes

-   -   Using probes of variable lengths opens the door for considerable        redundancy. Redundancy is desirable because it permits multiple        measurements to cross-check, but it may complicate matters by        altering the probability model surrounding a single Record.        Specifically, if a single Record represents two probe-landings        (one of “aaaaaaaa” and another of “aaaaaaaaa”) it should be        treated as though it has less error than would a single        probe-landing. To avoid these complications and keep the number        of reactions down, the following condition may be imposed on any        probe set: Let T denote the set of all infinitely-long strings        over the alphabet {“a,” “c,” “g,” “t”}. For all strings str in        T, there exists exactly one probe p of length l in our probe set        such that the l−1-length prefix o fp is a suffix of str. This is        identical to the condition that no probe in our probe set is the        suffix of another.

In certain embodiments, the disclosed methods for biopolymer sequencingare implemented as a computer program product for use with a computersystem. Such implementations may include a series of computerinstructions fixed either on a tangible medium, such as a computerreadable medium (e.g., a diskette, CD-ROM, ROM, or fixed disk) ortransmittable to a computer system, via a modem or other interfacedevice, such as a communications adapter connected to a network over amedium. The medium may be either a tangible medium (e.g., optical oranalog communications lines) or a medium implemented with wirelesstechniques (e.g., microwave, infrared or other transmission techniques).The series of computer instructions embodies all or part of thefunctionality previously described herein with respect to the system.Such computer instructions can be written in a number of programminglanguages for use with many computer architectures or operating systems.

Furthermore, such instructions may be stored in any memory device, suchas semiconductor, magnetic, optical or other memory devices, and may betransmitted using any communications technology, such as optical,infrared, microwave, or other transmission technologies. It is expectedthat such a computer program product may be distributed as a removablemedium with accompanying printed or electronic documentation (e.g.,shrink wrapped software), preloaded with a computer system (e.g., onsystem ROM or fixed disk), or distributed from a server or electronicbulletin board over the network (e.g., the Internet or World Wide Web).Some embodiments of the invention may be implemented as a combination ofboth software (e.g., a computer program product) and hardware. Stillother embodiments of the invention are implemented as entirely hardware,or entirely software (e.g., a computer program product).

In one embodiment, the steps described above are used to determine awhole or partial sequence of a target biopolymer. The method may includeproviding a double-stranded biopolymer target molecule. Then, the targetmolecule or fragment thereof may be contacted with a plurality of nprobes each having specificity for one or more recognition sites of thetarget molecule, thereby forming local ternary complexes along thedouble-stranded biopolymer target molecule or fragment thereof, each ofthe recognition sites having a known sequence. At least a subset of thelocal ternary complexes along the double-stranded biopolymer targetmolecule or fragment thereof may be detected to determine one or moreprobe maps. A whole or partial sequence of the target biopolymer may bedetermined using at least one or more probe maps.

Detection may include passing the double-stranded biopolymer targetmolecule or fragment thereof through a nanopore and detecting anelectrical signal indicative of locations of the local ternary complexesalong a length of the double-stranded biopolymer target molecule,thereby determining the one or more probe maps. This detection step maybe repeated for each of the n probes.

References, patents, and published patent applications specified hereinare incorporated herein by reference.

The described embodiments of the invention are intended to be merelyexemplary and numerous variations and modifications will be apparent tothose skilled in the art. All such variations and modifications areintended to be within the scope of the present invention as defined inthe appended claims.

What is claimed is:
 1. A method for determining a whole or partialsequence of a target biopolymer, the method comprising the steps of:providing a double-stranded biopolymer target molecule; contacting saidtarget molecule or fragment thereof with a first probe havingspecificity for one or more first recognition sites of said targetmolecule, thereby forming one or more first local ternary complexesalong said double-stranded biopolymer target molecule or fragmentthereof, said one or more first recognition sites each having a firstknown sequence, wherein the first probe is a member selected from thegroup consisting of an oligonucleotide, a peptide, and a polyamide;contacting said target molecule or fragment thereof with a second probehaving specificity for one or more second recognition sites of saidtarget molecule, thereby forming one or more second local ternarycomplexes along said double-stranded biopolymer target molecule orfragment thereof, said one or more second recognition sites each havinga second known sequence; passing said double-stranded biopolymer targetmolecule or fragment thereof through a detector configured to detect aternary complex through an increase or decrease of signal; detecting anelectrical signal indicative of one or more locations of at least asubset of said one or more first local ternary complexes along saiddouble-stranded biopolymer target molecule or fragment thereof todetermine a first probe map; detecting an electrical signal indicativeof one or more locations of at least a subset of said one or more secondlocal ternary complexes along said double-stranded biopolymer targetmolecule or fragment thereof to determine a second probe map; anddetermining said whole or partial sequence of said target biopolymerusing at least said first probe map and said second probe map.
 2. Amethod for determining a whole or partial sequence of a targetbiopolymer, the method comprising the steps of: providing asingle-stranded molecule; providing a double-stranded biopolymer targetmolecule by adding a second strand to the single-stranded molecule;contacting said target molecule or fragment thereof with a plurality ofn probes each having specificity for one or more recognition sites ofsaid target molecule, thereby forming local ternary complexes along saiddouble-stranded biopolymer target molecule or fragment thereof, each ofsaid recognition sites having a known sequence, wherein each of the nprobes is individually a member selected from the group consisting of anoligonucleotide, a peptide, and a polyamide; detecting at least a subsetof said local ternary complexes along said double-stranded biopolymertarget molecule or fragment thereof to determine one or more probe maps;and determining said whole or partial sequence of said target biopolymerusing at least said one or more probe maps.
 3. A method for determininga whole or partial sequence of a target biopolymer, the methodcomprising the steps of: providing a double-stranded biopolymer targetmolecule; contacting said target molecule or fragment thereof with aplurality of n probes each having specificity for one or morerecognition sites of said target molecule, thereby forming local ternarycomplexes along said double-stranded biopolymer target molecule orfragment thereof, each of said recognition sites having a knownsequence, wherein each of the n probes is individually a member selectedfrom the group consisting of an oligonucleotide, a peptide, and apolyamide; passing said double-stranded biopolymer target molecule orfragment thereof through a detector configured to detect a ternarycomplex through an increase or decrease of signal; detecting at least asubset of said local ternary complexes along said double-strandedbiopolymer target molecule or fragment thereof to determine one or moreprobe maps; and determining said whole or partial sequence of saidtarget biopolymer using at least said one or more probe maps.